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Abstract 

We present a perturbative QCD calculation of diffractive Dg production in 
charged current deep inelastic scattering. In the two-gluon exchange model, we 
analyze the diffractive process lyN —> which may provide useful informa¬ 

tion for the gluon structure of nucleons and the diffraction mechanism in QCD. The 
cross section of diffractive Dg production with xbj = 0.005-0.05 and Ei, = 50 Gev is 
found to be 2.7 x 10“® pb. In spite of this small cross section, the high luminosity 
available at the zz-Factory in the future would lead to a sizable number of diffraction 
events. 

PACS numbers: 42.25.Fx, 13.60.Le 

Diffractive leptoproduction of mesons has received much attention||l[|-P]due to two rea¬ 
sons. First it is of interest for the study of diffractive production mechanism within QCD 
and second, its cross section is dominantly proportional to the square of the gluon density 
in the nucleon, e.g., in the case of diffractive J/'ip electroproduction. 

Aside from the diffractive electroproduction processes, the charged current(CC) in¬ 
duced diffraction may also be interesting. To the lowest order in perturbative QCD, CC 
diffractive deep inelastic scattering (DIS) @] proceeds by the Cabibbo-favored production 
of the {ud) and (cs) states, and the two-gluon exchange between the (cs) and the nu¬ 
cleon may be the dominant mechanism for the diffractive production of charmed strange 
mesons. With the help of high luminosity available at the z/-Factory, neutrino induced 
diffraction in CC DIS can shed more light on the QCD mechanism of diffractive meson 
production. At the same time, it is a new way to study the gluon structure of nucleons. 

We now consider the diffractive process(Fig. |l|), 

+ V ^ /i- + V' + D+. (1) 

We shall be concerned with the kinematic region where Bjorken variable xbj = ■ q 

is small. The three-fold differential cross section is 

dcr _ _ 

dxBjdQ^dt 32(47r)3 sin^ 6w xs^s'^iQ'^ + 

where s = {p + lY k, 2p ■ I, t = vP and = —Q^, p and I are the 4-momentum of the 
nucleon and the lepton respectively. The leptonic tensor is 

L^,u = ( 3 ) 
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Figure 1: Diagram for the neutrino induced diffractive D+ production. 

To the lowest order in perturbative QCD, the hadronic current can be calculated 
from the colorless two-gluon exchange subprocesses shown in Fig. |^. We will use the 
nonrelativistic approximation writing the vertex in the form + M 25 J 75 . The 

constant specihes the cs coupling to the . We choose the wave function as 

kt) = {kt)6{z - ttIc/Md,), (4) 

where z and z = (1 — z) denote the fractions of Df momentum carried by the c and s 
quarks respectively, k is their relative momentum. Here we take k. irtc + mg- 

We hrst evaluate the gluon loops in the Feynman diagrams shown in Fig. ^ It is 
convenient to perform the loop integration in terms of Sudakov variables0. That is, for 
all particles the 4-momentum are decomposed in the form 

ki = aiq' + Pip' + kiT, (5) 

where p' and q' are respectively the light-like momenta of the nucleon and boson, 
that is, = 0. In particular 

p = p' + a^q', q = q' + (6) 

with = m?^/2p' ■ q' and = —Q‘^/2p' ■ q'. We consider the limit 2p' ■ q' S> 
then we have 2p' ■ q' k. 2p ■ q. 

Within the nonrelativistic approximation the quarks with momenta(see Fig. zq^ + 
K and zq^ — k are almost on mass shell. The integration over the gluon longitudinal 
momentum leads, in the hrst diagram, the upper quark with momentum zq^ — k + k to 
be on shell, leaving only the quark propagator {rf — to be integrated over in the 

gluon integration. 

Using the Sudakov decomposition, we hnd 

- ml = ~ (Q‘^ + k"^) (7) 

where Q = z{l — z){Q‘^ A which is the relevant effective perturbative QCD factor¬ 

ization scale [ 0 . 
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zq^-k+K zq^-\-K 



(a) 


zq^+K 



(b) 


Figure 2: Two of the four sub-diagrams. The other two diagrams are obtained by inter¬ 
changing c and s quark lines. 


Taking the CKM matrix element 14s = 1, we write the contribution given by the 
Feynman graph in Fig. |^a as 

/■d/cr^r[75(^^ + Mz3j^(^^^-^ + mc)7^(l-75)(/i+ms)^]0(fc) 


8v^2;7r sin 9 


w 


(2p ■ q)k‘^{k — uy{rf — m^) 


where Fc = 2/3 is the color factor, and (p{k) describes the emission of the gluon pair by 
the proton 0, 


Stt 

0 (^) = -^fBFKL{x,k‘^). 


(9) 


where is the gluon density unintegrated over k^ that satishes the BFKL equation 

which effectively resums the leading Q!s1ii[1/x] contributions, with 


X f3D, - f3^ = 


2p ■ q 


( 10 ) 


To relate fgp^^ to the conventional gluon density, which satishes GLAP evolution, we 
must integrate over k‘i 


xg{x, Q^) = J {x, kl) 


( 11 ) 


In analogy to the derivation of (^, we hnd the sum of the four diagrams in Fig. ^ is 




3ir!z(l - j)eg„g;n(2p ■ g) r dk-p q d(xg(x,k^)) 


i=i sin 9w 

To the lowest order in k^, we have 

?>mz{l - z)eggglFc{2p ■ q)q^ 


Q\q' + dkl 


( 12 ) 




[xgix, Q^)). 


^4 


3/2 sin 9wQ 


(13) 


So far we have calculated only the imaginary part of the amplitude. We can use 
dispersion relations 0, O to determine the real part, and numerically we hnd it to be 
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Figure 3: Differential cross sections as a function of xbj and Q^. The neutrino energy has 
been chosen as = 50 Gev (solid lines) and Ei, = 40 Gev (dotted lines). 


not negligible. Including the real part contribution as a perturbation we now rewrite the 
differential cross section (^) as 


dcr e'^{s — 2p ■ q)\A\‘^ 

dxBjdQ^dt 16(47r)3 sin^ 


(14) 


where 


127 r‘^ieg^as{Q^)Ec. ^^2. , in d{xg{x,Q'^)) 

= - - -=0— \xq{x,Q ) H-—- 

■^siaewQ 2 aim 


(15) 


Since we are concerned with small x, the effect of the nonzero value of \t\, of which 
the minimum is is expected to be small. Then we can integrate out t by 


ydte-“ = ^ (16) 

where we will use the experimental slope value b = 3.3 Gev“^ as in similar processes]^. 

To give numerical results, we take the input parameters as follows: Mw = 80.4 Gev, 
rric = 1.5 Gev, rUg = 0.5 Gev. The running strong coupling constant is chosen with 
as(m^) = 0.27. For the gluon distribution function, we select the Gluck-Reya-Vogt(GRV) 
next-to-leading order(NLO) set|@]. The constant 51 ^ can be expressed in terms of the 
decay constant by 


(0|c7/,(l - 75)s|^^) = 

which gives g^ = /d„/ 4. Here we choose fos = 280 Mev^]. 


( 17 ) 
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In Fig. ^(solid lines) we show the resnlts obtained for the differential cross sections 
dcr/d^Bj and dcr/dQ^. The nentrino energy has been chosen as E^, = 50 Gev. For the 
plot of XBj-dependence, has been integrated from 0.5 Gev^ to the upper bound given 
by the constraint on inelasticity y = Q^/2a:BjP ■ / < 1. In the plot of Q^-dependence, xbj 
has been integrated from lower bound to 0.05 and taking the same kinematic constraint 
mentioned above. Integrating over and xbj in the kinematical region specihed above 
gives a value for the total cross section of a = 2.7 x 10-5 pb. 

To see the sensitivity of the differential cross sections to the neutrino energy, we also 
present the results for E^, = 40 Gev in Fig. |]( dotted lines). The kinematic regions of 
and XBj are the same as in the Ey = 50 Gev case except that the upper bound of is 
chosen as 0.065. Integrating out all variables gives the total cross section of cr = 2.0 x lO-^ 
pb. In spite of the small cross section, the high luminosity available at the n-factory in 
the future would lead to a sizable number of events of the order of magnitude 10^. 

Some discussions are in order. First, in the two gluons exchange processes in general 
we should encounter the so-called off-diagonal gluon distribution function |^. But it is 
expected that for small x there is no big difference between the off-diagonal and the usual 
diagonal gluon densities 0- So in the above calculations we have estimated the small 
X production rate by approximating the off-diagonal gluon density by the usual gluon 
density. This situation is similar to many diffractive production processes at hadron 
colliders |T3|, |T4|(for detailed discussions, see IT^). 

Second, we have used Q = z{l — z){Q‘^ -|-as the energy scale for the application 
of the perturbative QGD. The applicability of pQGD is guaranteed by the large value of 
^ {rric -l- rris)'^. So can be chosen to be rather small, say, 0.5 ~ 1.0 Gev^. 

Third, we have used nonrelativistic approximation to describe the Dg wavefunction, 
and this will cause some uncertainties in our calculation. Relativistic effects can be quite 
important and should be further considered in a similar way as in Big. 

In conclusion, we have calculated the diffractive Dg production rate in the neutrino 
induced charge current DIS process in the two-gluon exchange model in QGD, and found 
the diffractive production of Dg to be observable with the high luminosity available at the 
i/-Factory in the future. 

After the calculation in this work was completed, B. Lehmann-Dronke and A. Schafer 
[T^ published a preprint treating a similar process to that we considered. But they 


analyzed exclusive Dg production in the large xbj region, whereas we studied the diffractive 
production of the Dg meson in the two-gluon exchange model with small xbj- Although 
calculated in different methods and in different kinematic regions, our total cross section 
has the same magnitude as theirs. 

We would like to thank F. Yuan for his valuable discussions and K.Y. Liu for some 
numerical calculations. This work was supported in part by the National Natural Science 
Foundation of P.R. Ghina, and the Education Ministry of P.R. Ghina. 
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